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ANDRE REZNIKOV 

Abstract. We consider the hyperbolic Casimir operator C defined on the 
tangent sphere bundle SY of a compact hyperbolic Riemann surface Y. We 
prove a nontrivial bound on the L 2 -norm of the restriction of eigenfunctions 
of C to certain natural hypersurfaces in SY. The result that we obtain goes 
beyond known (sharp) local bounds of L. Hormander. 



1. Introduction 

1.1. Restriction problem. The purpose of this paper is to draw attention to the 
following well-known question in PDE. Let X be a compact n-dimensional smooth 
manifold, and let C be a second order hyperbolic operator (i.e., an operator with 
the principal symbol p of the signature (n — 1, 1)). Let O C X be a closed smooth 
hypersurface, which we assume to be non-characteristic for C . Let / £ C°°(X) be 
a smooth function. What can be said about the norm of the restriction of / to O 
in terms of norms of / and C(f) on XI 

The "local" theory related to such questions goes back to works of L. Hormander 
(for example, see [HI Theorem 25.3.11]). One knows that the answer depends on 
the curvature of O with respect to the pseudometric defined by the principal 
symbol p. Resulting bounds and various extensions are discussed at length in the 
paper by D. Tataru, [Taj . In particular, for the smooth time-like hypersurface we 
have 

(1-1) ll/|o|| w 3/. (0) < A ■ (||/||v^ c (x) + \\C(f)\\ Llc{x) ) , 

for some constant A > depending on the geometry of X, O and p. Here we 
denote by XV s the L 2 -Sobolev norm of the order s on X (for example, associated 
to a choice of Riemannian metric on X) . This should be compared with the Sobolev 
embedding theorem giving the loss of 1/2 of the derivative (i.e., ||/|o|| T/I/ i/2/ n >, < 

w ioc (y) 

A' ■ WfWw, 1 (x))- If m addition O is curved, A. Greenleaf and A. Seeger [GrSj 

obtained the improvement ll/|o|| w s/6 (0) < A" ■ (\\f\\wt oe (x) + W C (f)\\Lf oc (X)) of 
the bound (jl.ip . These bounds are optimal, as could be seen from the constant 
coefficients examples. 

Since these results are local, one can not hope to improve these bounds for gen- 
eral functions using global geometry of X, O and p. However, we want to present 
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an example where the improvement is possible for functions having dominating 
first term on the right in the bound (|1.1[) . Namely, we will construct an example 
where 

(1-2) \\f\o\\ L -{o)<A a>e -\\f\\ W s {x) , 

for any e > 0, and for any / £ C°°{X) satisfying ||C(/)||iy e (x) < a for a given 
a > 0. In particular, for eigenfunctions of C with bounded eigenvalue, there 
is essentially no loss of smoothness in taking the restriction (at least to some 
special hypersurfaces) . This raises the question if such a phenomenon persists 
more generally. Namely, one would like to see if it is possible to impose certain 
conditions on the geometry of X, O and p which would imply that 

(1.3) \\f\o\\mo) < A x ,o, Pt e ■ (WfWw'W + \\C(f)\\wHx)) 

for eigenfunctions of C. 

As we mentioned above, one can not hope for a local explanation for a bound 
of the type (|1.3I) . In fact, one can easily construct an example (based on spherical 
harmonics, as usual) where nothing close to the bound (jl.3[) holds for a compact 
X and O (see Example 11.51 below). Nevertheless, an intuition from recent works 
on Quantum Chaos suggests that there are some natural situations where one 
might expect some version of (II .3|) to hold. This paper proposes one such an 
example which is based on the theory of hyperbolic Riemann surfaces (or more 
appropriately, on representation theory of PGL2(M)). We show that in that par- 
ticular example the bound (11.31) holds for the "low energy" spectrum, and for 
special hypersurfaces O (e.g., associated to closed geodesic circles on the Riemann 
surface). 

Our proof (and the example) comes from representation theory, and is a mixture 
of results and techniques from |Re2) and |Re3j . 

1.2. Casimir operator. Here we describe our construction. Namely, we con- 
struct a 3-dimensional compact manifold X (S^-fibered over a compact hyperbolic 
Riemann surface Y), and a nondegenerate second order hyperbolic operator C on 
X (the Casimir operator) such that for any immersed two-dimensional tori O d X 
which is fibered over a closed geodesic or a closed geodesic circle on Y , the bound 
(0 holds. 

First, let us recall the well-known setup (due to I. Gelfand and S. Fomin, |G6j ) 
which allows one to apply the representation theory to some special PDE operators. 

1.2.1. Hyperbolic geometry: Riemann surfaces and Geodesies. Consider a compact 
Riemann surface (Y,^) endowed with the constant negative curvature metric 172, 
the corresponding volume element dv2, and let A be the corresponding (nonneg- 
ative) Laplace operator on Y. Let SY be the tangent sphere bundle over Y. 
It is more convenient for our purposes to work with the oriented orthonormal 
frame bundle RY over Y. The 3-dimensional manifold RY has two connected 
components for the oriented Riemann surface Y, and there is the projection map 
RY — > SY taking the first vector in the frame. 

The central object in our method is the group G = PGL2 (R) (one can work with 
the group SL^K) of matrices g with det(g) = ±1 instead). This group naturally 
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appears in geometry as the group of isometries of the hyperbolic plane H not 
necessarily preserving orientation. We fix the standard maximal compact subgroup 
K = PO(2) C G, and denote by K = PSO(2) C PGL^(K) its maximal connected 
subgroup (Kg ~ S 1 ). We have the isomorphism T-L ~ G/K. The uniformization 
theorem implies that there exists a lattice L C G (i.e., a discrete co-compact 
subgroup r ~ 7Ti(F)) such that Y = T\H. We set X = T \ G. This will be our 
manifold. It is important to note that X is a homogeneous manifold for G acting on 
the right. We have the isomorphism X ~ RY. We will assume for simplicity that 
r C PSL2(R) (i.e., we distinguish between Y and its complex conjugate Riemann 
surface Y). In that case, X consists of two connected components. While it seems 
artificial at first, the disconnectedness of G is important in our considerations (and 
this is the reason we work with PGL2(R) and not with SLa(K)). 

The Ricmannian metric <? 2 could be extended to the Riemannian metric 173 on 
X. In fact, it is easier to choose (73 first by choosing an Ad(i^o)-i nv ariant form on 
the Lie algebra q = s/2 (which is isomorphic to the tangent space at any point of 
X), and transporting it to the whole X by the action of G. We also fix the (unique 
up to a constant) G-invariant measure dx on X (e.g., normalized by Vol(X) = 1). 
We denote by p : X —> Y the corresponding projection and by B x ~ S 1 x Z 2 the 
fiber passing through x £ X (in fact, B x = xK is an Zf-orbit). 

Consider a closed geodesic I C Y with the natural line element dl on it, and 
the 2-dimensional tori O = p~ l {l) C X (O ~ I x S 1 ) obtained as the preimage of 
I under p (in general this will consist of two copies of 2-tori). Let T C G be the 
subgroup of diagonal matrices. It is well known that I is an orbit of T (in fact, 
any closed orbit of T gives rise to a closed geodesic on Y) . Hence, O = I ■ K c X 
is the result of the action of the compact subgroup K on the set I. Integrating dl 
with respect to this action, we obtain an area element do on O. 

The 3-dimensional manifold X and the hypersurface O C X described above 
are the geometric ingredients of our example. In fact, we will change the setup in 
the proof in order to avoid certain technical complications explained in Remark 
12. 5. II Instead we will deal with another family of tori in X. These will also 
be S' 1 -fibered, but this time over a geodesic circle on Y. The case of a closed 
geodesic will be discussed elsewhere. Hence, in this paper, we present the proof 
for a tori O = p^ 1 ^) C X fibered over a geodesic circle a C Y, a — cr yo (r) — {y € 
Y I d(y,yo) — r} centered at an arbitrary point j/o G Y. In fact, we will choose a 
connected component of this set. The resulting (immersed) hypersurface O clearly 
is _ftTo-invariant (topologically it is a 2-dimensional tori). 

1.2.2. Casimir operator. There exists a unique (up to a multiplicative constant) 
second order differential operator on X which is invariant under the (right) action 
of G on X. We choose such an operator G, called the Casimir operator on X. In 
physics it naturally appears in the Kaluza-Klein theory (see [Z]). 

Let db — d/db be the unit vector field on X tangential to fibers of the projection 
p : X — > Y. We choose G (or A) by requiring that G = A + &l (or rather 
requiring that the G-invariant operator G coincides with A on functions constant 

along fibers of p). The symbol of G has the signature (+ H ). We also will 

need an elliptic operator on X in order to define Sobolcv norms. We will use 
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Ax = A — d 2 = C — 2d 2 which coincides (under the appropriate normalization) 
with the Laplacian on X induced by the Riemannian metric (73. 

1.2.3. Eigenf unctions. The operator C on X has a nice spectral theory (which is 
admittedly unusual for a hyperbolic operator). Namely, the spectrum of C coin- 
cides with that of the elliptic operator A (allowed also to act on forms). However, 
while all eigenspaces of A arc hnitc-dimcnsional, all nontrivial eigenspaces of C 
are infinite-dimensional. 

One can interpret eigenfunctions of C in terms of the representation theory of 
PGL2(M). We denote by R(g)f(x) = f(xg) the action of G on functions on X. 
Let A £ R be an eigenvalue of C. In what follows, it will be more natural (from 
the point of view of representation theory) to write A = (1 — A 2 )/4 for appropriate 
A G C. Let U\ C C°°(X) be the space of smooth eigenfunctions of C with the 
eigenvalue A. 

1.3. The restriction. We fix the set O C X as above (i.e., we fix a point yo G Y, 
and the radius r > of the geodesic circle <r VQ (r); O = P~ 1 {o~ Vo (r))). We are 
interested in the norm of the natural restriction map on the subspace Ua- Namely, 
consider the restriction map ro,A ■ Ua ~ > C°°(0), given by ro,A(/) = f\o f°r 
/ G Ua, and the corresponding Hermitian form 

(1-4) H A {f) = [ \roAf)\ 2 do 

Jo 

defined on the space Ua- The form Ha is ifo-invariant: H\(R(k)f) = H\(f) for 
all k G Kq, since the set O and the measure do are ifo-iuvariant. 

Our main result is 

1.4. Theorem. For any e > 0, there exists a constant Aa. £ > such that 
(1-5) H A (f)<A A , e -\\f\\ 2 Wsix) , 

for any f G Ua- 

Here we denote by W S (X) the L 2 -Sobolev norm of order s on X (see |BRlj for 
the representation theoretic treatment). To define these norms, one can use the 
elliptic operator A x on X (e.g., ||/||^ a(x) = \\f\\h {x) + ||A X (/)||| 2(X) ). The 
value s = corresponds to the L 2 -norm. 

Our approach to Theorem 11.41 is based on the celebrated theorem of Gelfand 
and Fomin |G6j . which interprets spaces Ua in terms of irreducible unitary rep- 
resentations of PGL2(M) appearing in L 2 (X,dx) (i.e., each space Ua is a finite 
direct sum of isomorphic irreducible representations of G; see Section [5]). These 
irreducible representations are called automorphic representations. This allows 
one to use powerful methods of the representation theory. We deduce the bound 
(jl.5p from the corresponding bound for each irreducible component inside of the 
space Ua (see Theorem 12 .7[) . 

Remarks. 1. The bound (|1.5[) means that on an eigenspace Ua the restriction 
map ro.A is L 2 — W £ bounded for any e > 0. We note, however, that the bound 
we obtain for the constant Aa. 6 (i.e., the norm of the restriction map) is quite 
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poor, and corresponds to the general bound for the norm of the restriction of 
eigcnfunctions to geodesies on Y (i.e., it does not distinguish between positive and 
negative curvature). Quantum Chaos intuition suggests the following 

Conjecture. For any fixed e > 0, 

A A , e « |A| £ 

for any e > 0. 

This conjecture is equivalent to the bound (jl.3[) for our particular X and O. 

2. The space Ua has the natural structure of the tensor product V\®M\, where 
V\ is an irreducible (infinite-dimensional for A / 0) unitarizable representation 
of G, and M\ is a finite dimensional complex vector space with the dimension 
dim M\ = m\ equal to the multiplicity of the eigenvalue A of A on Y (one can 
take M\ to be the A-eigenspace of A on 7). We deal with the restriction norm 
for a single space V\, and them use the fact that M\ is finite-dimensional. This is 
equivalent to choosing an orthonormal basis in the space M\, or what is the same, 
choosing an orthonormal basis in the A-eigenspace of A. 

3. A more interesting form to consider is the form which is obtained as the 
restriction to a single geodesic l; namely, the Hermitian form h\(f) = J, \f\i\ 2 dl 
on U\. The form h\ is not L 2 -bounded. Using a pointwise bound (e.g., as in 
[BR1 ), it is easy to see that it satisfy the bound h\(f) < A e \\f\\'^ vl/2+E ^ x ^ for 
any e > and some A e . Forms H\ and h\ are related via integration over Kq: 
H\ = J Ka k ■ h^dk, and hence it is not surprising that the form Ha is "smoother" . 

4. It seems that our method should provide the following information about 
the constant Aa, e (at least for O fibered over a geodesic circle a). 

Let Ab = —dgi be the Laplacian along S 1 fibers of the bundle O — > a. The 
operator Ab commutes with C since it corresponds to the action of PSO(2) on X, 
and has the simple spectrum n 2 , n € 2Z. We denote by H\ the form H\ restricted 
to the (irreducible) automorphic representation V\. Hence for every eigenspace V\, 
we can choose an orthonormal with respect to the form H\ basis {</>™} n e2Z of joint 
eigenfunctions of C and Ab (C<^ = A<^, Ab4>^ = n 2 </)"). Functions cj) 7 ^ are 
exponents along the fibers S . Such a decomposition follows from representation 
theory. 

To bound the norm of the form H\, one needs to bound coefficients Hx^ 1 ^). 
It is plausible that with more work the following bound could be extracted from 
our method: -Ha(</>™) — C| n | e > for \n\ > |A| and some universal C > 0. This would 
mean that there exists a universal constant in (|1.5|) for \n\ > |A|. For \n\ < |A|, 
we have essentially the standard bound of Hormander (or rather the improvement 
for the restriction to a curved hypersurface) H\((f?^) < IAI 1 / 6 , i.e., A e ^\ < | | x / e in 
this range. 

One might expect that orthogonal irreducible representations inside the space 
Ua are "nearly" orthogonal with respect to the form Ha, but we do not know 
how to prove it. In fact the multiplicities a less problematic for our method than 
"clusters" of eigenvalues (and for eigenpackets one expects a different answer). 

As a side remark, we note that there is a natural setup of "ladders" which allows 
one to consider eigcnfunctions 0" in the framework of PDE theory (in fact, the 
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notion exists on any Riemannian manifold). In particular, Guillemin, Sternberg 
and Uribe introduced a notion of "fuzzy ladders" (see |GSj . |GUj . |ST) . [Z])- In our 
setting, this notion concerns a sequence of joint eigenfunctions of two commuting 
operators: 8b (the Kaluza-Klein operator) and Ax on X. One sees immediately 
that these joint eigenfunctions are nothing else but the eigenfunctions (/>". 

We note that there is a variety of effective results about the asymptotic behavior 
of eigenfunctions (f>^ (e.g., effective equidistribution of |<^™| 2 on SY as |n| — > oo 
shown in [Rel , and uniform bounds for ||<^||l 4 (x) for |n| > |A| mentioned in 
[BR3 ). These results are obtained through representation theory methods. 

5. Results discussed in this paper could be easily extended to noncompact 
hyperbolic surfaces of finite volume. In fact, this what generated our original 
interest in the problem. For such surfaces of the arithmetic origin (e.g., Y = 
PSL 2 (Z) \ W), the Theorem 11.41 is equivalent to a subconvexity bound on certain 
L-functions from Number Theory (sec Re3 ). In particular, the corresponding 
bound should hold for the Eisenstein series, and also for horocycles. 

1.5. Examples related to S* 2 . For some pairs O C X, the bound (jl.3p does not 
hold. We briefly discuss some examples related to S 1 bundles over the standard 
sphere S 2 . 

First we give a simple example where the local bound (jl.ip is sharp. Consider 
X = S 2 x S\ O = (equator of S 2 ) x S 1 C X, and C = A S 2 - |A S i (i.e., X 
could be viewed as a wrapped manifold; here As» is the (positive) Laplacian on 
the standard sphere). The hypersurface O is non-characteristic, time-like and 
flat. An easy calculation shows that the eigenspace V1/4 of C with the eigenvalue 
A = 1/4 is infinite-dimensional. An orthonormal basis of V1/4 consists of functions 
y m ,i{ s it) — c-Y^ m \s)e^ 2l+lS>t , s £ S 2 andi € S 1 . Here are norm one classical 
spherical harmonics with the eigenvalue 1(1 + 1) (i.e., A32Y, ' — 1(1 + l)Y^ m \ 
\m\ < I), and c = l/y/2n is a normalizing constant. It is easy to see that the local 
bound p. ID is achieved in this example for functions in the infinite-dimensional 
space V1/4. This is due to the fact that there are spherical harmonics having 
maximal possible L 2 -norm on the equator. 

A more relevant example to our discussion would be the tangent sphere bundle 
X = S(S 2 ) over S 2 . Here the situation is more complicated. Again we can 
use representation theory since X ~ SO(3). Hence we choose X = SO(3) and 
G = SO (3) acting on X on the right (we also use the left action below). Let 
/ C S 2 be a circle on S 2 and O its preimage in X (O is a 2-dimensional tori as 
before). We can view I as an orbit of a compact subgroup K' ~ S 1 C SO (3) acting 
on the left on X. Hence O is a double K' x SO(2) orbit in X (in such a picture 
K' could coincide with SO(2)). 

There is no natural hyperbolic operator on X from the point of view of repre- 
sentation theory, but there is a family of such operators which could be described 
via the right action of SO (2) on X. Let C3 be the (elliptic) Casimir operator on 
X (i.e., the unique (up to a multiple) second order SO(3)-invariant differential 
operator on X; we normalize it to coincide with the Laplacian A2 on 5* 2 ), and 
let C2 be the differential operator associated with the Casimir operator of SO (2) 
(i.e., coming from the action of SO(2) on X). We consider a family of operators 
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Ct = C3 — t ■ C<x depending on the real parameter tel. To obtain a hyperbolic 
operator, we have to choose t > 1. The distribution of eigenvalues of C t depends 
on the size of t and on its Diophantine property. In particular, for t which is 
irrational, one can see that for big enough T, the space Vt which is generated by 
eigenfunctions of C t with the absolute value of the eigenvalue bounded by T is 
infinite-dimensional. The same is true for a big enough rational t. We can use the 
Peter- Weyl theorem on the structure of L 2 (SO(3)) in order to describe the space 
Vt- It is not difficult to see that Vt is the direct sum of irreducible representations 
of SO(3), and dimensions of these representations are unbounded. The study of 
the corresponding restriction problem could be reduced to asymptotics of gener- 
alized spherical functions. One can see then that there is an infinite sequence of 
eigenfunctions in Vt having on O polynomially big (with respect to their Sobolev 
norm) restriction norm. However, the exact rate of this growth and its dependence 
on O and t is more difficult to determine and deserves further investigation. 

Finally, we note that the bound (|I.3[) holds for a flat tori T 2 C T 3 and the 
appropriate (constant coefficients) hyperbolic operator. 



1.6. Idea of the proof of Theorem 11.41 The proof of bound (jl.5D is based on 
techniques from representation theory. However, the basic idea is quite elementary. 
We discuss it for a closed geodesic I, while in practice we give the proof for a 
geodesic circle and will discuss closed geodesies elsewhere. 

Let I C O C X be a closed geodesic. In particular, I = x$T c X is a closed 
orbit for the diagonal subgroup T. Let V\ C U\ be an irreducible representation 
of G. We denote by H\ the form H\ restricted to the automorphic irreducible 
representation V\. We bound the form H\, i.e., the form restricted to any 
irreducible component inside of U\. The space U\ is isomorphic, as an abstract 
representation of G, to a finite number of m\ copies of V\. To obtain the bound 
for the form H\, we use the fact that m\ is finite. A simple reduction shows 
(see section \2. 61) that in order to bound the norm of the Hermitian form H\, it is 
enough to bound values of the form hx(4>e n ) = Jr I0e„| 2 ^ for some special vectors 
4>e n € V\ (these are vectors which are exponents e mB along fibers B x ~ S 1 ; see 
section 1272]) . Here we use the crucial fact that H\ (and H\) is i^-invariant, and 
the structure of irreducible representations of G (i.e., that the space of so-called 
-fTcrtypes are one-dimensional). 

Hence we are interested in values of h\(<fi en ) — f, \<j) en \ 2 dl for |n| — > 00. We 
have I ~ S 1 (in general I ~ S 1 U S 1 is a union of two closed geodesies, but we 
disregard this complication since in this paper we will only deal with geodesic 
circles where this complication could be avoided) . Hence we can use Fourier series 
expansion on I. Consider expansions <j> en \i(9) = 'Yl k ak{4 > e rl )e lke and h\(cf) en ) = 
S/c l a fc(^e„)| 2 - It turns out that representation theory implies that the coefficients 
a>k{4>e n ) could be naturally represented in the form ak{<fi e „) = a>k ■ h(k,ri), where 
the function h(k,n) is some special function which is well approximated by the 
classical Airy function. In particular, from such a representation it follows that 
for |fc| < (l-e)H, we have \h(k,n)\ 2 = 0(|ti| _:l ), and that for |fc| > (l + s)\n\, we 
have \h(k, n)| 2 = 0(\n\~ N ) for any N > and any e > 0. In the resonance regime, 
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k ~ n, we have the Airy type behavior \h(k, n)\ 2 ~ |n|~ 2 / 3 . We note that the 
coefficients cifc are well known in the theory of automorphic functions (see |Re3j ) . 

Hence we need to estimate the coefficients |afe| 2 . From simple geometric consid- 
erations, one can see that J2\k\<T \ a k\ 2 < A\T for any T > and some constant 
A\ > depending on A. This immediately implies that in the sum h\((f> en ) = 
12k l a fc(^e„)| 2 = J2k \ a k\ 2 \h(k, n )| 2 ? the contribution from regular ranges |fc| < 
(1 — e)\n\ and |fc| > (1 + s)\n\ is uniformly bounded as \n\ — > oo. 

Our main point then is that the bound 

(1.6) i flfe i 2 ^ B ^ Tl+£ 

\\k\-T\<ri 

holds for any T > and e > 0, and some constant B\ e > depending on A and 
e. Such a bound clearly implies our main claim (|1.5p in view of asymptotic for 
\h(k,n)\. To prove fj 1 . 6f) . we use another expansion of h\ and the related (less 
known) representation theory. 

Consider the following collection of sets: Al C I x 1 c I x I and Al C 
AX a X x X . The Hermitian form h\ could be viewed as a linear form on 
Va ® Vx C C°°(X x X) given by the integral over Al C X x X. Clearly, we 
used the action of T x T on I x i C X x X in order to describe the Fourier 
decomposition h\{4> en ) — J, |</> e J 2 cM = |afc((/) e?i )| 2 . We can use now the action 
of AG on AX C X x X in order to give another decomposition of the same form, 
hence leading to a certain useful identity. 

The second decomposition is also easy to describe in elementary terms as follows. 
Consider the function |</> e „| 2 € C°°(X). We have the spectral decomposition with 
respect to the Casimir operator: |0 e „| 2 = J2xieSpec(c) P r i(l^e„ | 2 ), where pn is 
the projection to the eigenspace V\ t . This implies the spectral decomposition 
h\(<t>e„) — Ji \4>e„ \ 2 dl — Y^,x JiP r i(\4>e n \ 2 )dl. It turns out that again representation 
theory allows one to identify quantities /, pri(\<f> en \ 2 )dl in some abstract terms (see 
Section[3]), and in particular again allows one to write J, pri(\<p en \ 2 )dl — b\ i -g{\i, n) 
as a product of an "arithmetic" coefficient bxi and a special function g(Xi, n). This 
shows that we have the identity 

(1.7) ^|a fe | 2 |/i(fe,n)| 2 = /i A (</» e J= J2 bxt-giKn). 

k X % eSpec(C) 

We use this identity to deduce the bound T | <T § \ a k\ 2 < B\^ £ Ti +e . 

The coefficients satisfy the mean- value bound Yl\x-\<T I^J 2 — CT 2 (which 
is in disguise a bound of L. Hormander for the spectral average for a value at a 
point of eigenfunctions of the Laplacian on Y). The function g(A, n) is again some 
kind of a special function with the Airy type behavior for |A| ~ n, and typical 
value of order of |A| _1 for non-resonance regime. It seems we are back to the same 
problem and have gained nothing. 

However, and this is the main reason we are able to show some non-trivial 
saving, the contribution to the Airy behavior of h(k, n) and that of g(X, n) comes 
from different sets of the "microlocalization" of the function (f> en (the meaning 
of this in representation theoretic terms is explained in section 12. 3[) . In other 
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words, the wavefront sets of distributions defining densities h(k, n) and g(X,n) on 
the space V\ are disjoint. This allows us to construct an appropriate test vector 
v £ V\ such that, on the one hand, v picks up the sum Y). nl 2 |afc| 2 on the 

left in the identity (|1.7p . and, on the other hand, has smooth behavior on the right 
part of (ll.7[) . The mean- value bound on coefficients b\ i then implies the bound 
CEU). 

2. Representation theory 

2.1 . Gelfand pairs. We will base our analysis on the notion of Gelfand pairs (see 
[GPJ. Instead of giving the general definition, we list two cases we need in the 
proof. 

• For any smooth irreducible representation V of G = PGL^R), the space 
Horrift- (V, y) is at most one-dimensional for any character x '■ Kq — > C x 
of a maximal connected compact subgroup of G (e.g., Kq = PSO(2)). 

• For any three smooth irreducible representations Vt, V2, V3 of G, the 
space of trilinear invariant functionals Homc(Vi (g> V2 <8> V3, C) is at most 
one-dimensional. 

In fact, we need only to consider unitarizable representations. We note that the 
second multiplicity one statement does not hold for SLi2(R) or PGL + (R) (the 
space is at most two dimensional). This is the reason we have to deal with the 
disconnectedness of PGL2(R). 

There is another well-known multiplicity one statement. For any smooth irre- 
ducible representation V of G, the space Homx(V,x) ls a t most one-dimensional 
for any character \ : T — > C x of a maximal tori in G (e.g., for the full diagonal 
subgroup in PGL 2 (M)). This is relevant to the discussion of closed geodesies. 

In order to set notations, we quickly review the standard constructions from 
the theory of automorphic functions and the relevant representation theory (see 
[G6], [Bo], pS], [La]). 

2.2. Irreducible representations of PGL^R). Infinite-dimensional irreducible 
unitary representations of G are naturally split into two types: induced representa- 
tions and discrete series. Unitary induced representations 7T\^ e are parameterized 
by a pair (A, e), where A is a complex number which belongs to the set £ «RU (0, 1) 
and e £ {0,1}. Discrete series representations tt^ are parameterized by an even 
positive integer k £ 2Z + . The parameter A we used before describes the action of 
a particular generator of the center of the universal enveloping algebra U(sl2) (i.e., 
the Casimir operator). The relation between two parameters is A = (1 — A 2 )/4 for 
induced representations, and A = fc/2(l — fc/2) for the integer case. Induced rep- 
resentations parameterized by iR are called principal series representations, and 
those parameterized by (0, 1) are called complementary series representations. We 
call these even for e = 0, and odd for e = 1 (although this deviates from the stan- 
dard terminology where this term is used in connection with the central character 
of GLi2(R)). Unitary induced representations and discrete series representations, 
together with two one-dimensional representations, det e , e £ {0, 1}, exhaust all 
irreducible unitary representations of G. 
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We note that only representations of principal and complementary series repre- 
sentations have a non-zero ATo-fixed vector, and hence give rise to eigenfunctions 
of A on Y. Representations of the complementary series correspond to what is 
called the exceptional spectrum of Y (i.e. A < 1/4). Representations of the dis- 
crete series correspond to forms on Y. We will assume for simplicity that there is 
no exceptional spectrum, i.e., that all automorphic representations with non-zero 
Ao-fixed vectors are representations of the unitary principal series. 

The structure of unitary representation of PGL2 (K) is well known. Induced rep- 
resentations could be modeled in various spaces of functions on various manifolds. 
In particular, we have the following realization of the space V\_ £ of smooth vectors 
for the induced representation 7Ta, £ . Consider the space Hx, s of smooth even homo- 
geneous functions on R 2 \0 of the homogeneous degree A— 1 (i.e., f(tv) = |t| A-1 /(i;) 
for any i G K x and ^ » £ I 2 ). The evenness condition is necessary since we 
consider the group PGL2(M). We have the natural action of GL^M) given by 
^X,e(g)f(v) = f{g~ x v) ■ I det(g)|( A_1 )/ 2 • det(g) £ , which is trivial on the center and 
hence defines a representation of PGL2(K). We call such a realization the plane re- 
alization. We will use the realization of such a representation in the space C^(S 1 ) 
of smooth even functions on S 1 (i.e., f(0 + ir) = /(#)) which is obtained by taking 
the restriction to the circle S 1 C R 2 \ 0. We call such a realization the circle 
model. In such a realization, the action of Kg ~ S 1 is given by the rotation. The 
invariant unitary norm then coincides with the standard norm on L 2 (S' 1 ) (this is 
where the assumption A G iR is used; for complimentary series the norm is not 
local). Hence there is a natural orthonormal basis {e n (9) = e m0 } n£ 2Z consisting 
of AT -equivariant vectors (called Ao-types). The parameter e specifies the action 
of the element 6 = ( ~~ x ) with the negative determinant. Namely, the action of 
S is given by nx, s {S)f(6) = (-l) e /(-6»). In particular, 7r A , £ (<5)e„ = (-l) £ e_„ for 
the basis of AT -types. 

The discrete series representations lack such a simple geometric model (a source 
of many computational difficulties) , and we discuss necessary amendments inlAl 



2.3. Automorphic representations. The notion of automorphic representation 
allows one to use effectively models of representations in our setup. Namely, till 
now we viewed the space V\ C C°°(X) as a subspace of functions on X. This is 
a rather inaccessible (however remarkable) realization of an irreducible represen- 
tation. For principal series, the alluded Gelfand and Fomin theorem implies that 
there is a G-equivariant map v\ : C^(S 1 ) —> V\, which we can assume to be an 
isometry. This means that every vector (ft G V\ is of the form <f> u = v\{u) for some 
function u £ C^(S 1 ). This allows us (sometimes) to translate various questions 
about eigenfunctions of C on X (or of A on Y) into questions about v\ and an 
appropriate vector u £ C^(S 1 ). We may argue that the correspondence cf> -f> u is 
somewhat analogous to what is called microlocalization of eigenfunctions. 

For discrete series, similar considerations apply. However, discrete series rep- 
resentations lack a good realization similar to the circle model above. We will 
discuss the appropriate changes in El 
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We assumed that the space X is compact (i..e, the Riemann surface Y is com- 
pact). Let 

(2-1) L 2 (X) = (®i(L Titet ,uj) ® (e^Lfe,,^)) 

be the decomposition into the orthogonal sum of irreducible unitary representa- 
tions of G. Here Vi : Li ~ L Tii£i — > £ 2 (A) are unitary representations of class one 
(i.e., those which correspond to Maass forms on Y with the eigenvalue (1 — rf )/4), 
and : Lfc — > L 2 (X) are representations of discrete series (i.e., those which corre- 
spond to holomorphic forms on Y). We denote by Vi C Li the corresponding spaces 
of smooth vectors and by pr £ . : L 2 (X) — > Vi(Li) the corresponding orthogonal pro- 
jections (note that pr £ . : C°°(X) — > Vi). Note that there might be multiplicities 
for automorphic representations (i.e., two realizations v, v' : — > L 2 (X) of an 
irreducible representation (L v ,n), having different images) notably for represen- 
tations of discrete series. 

2.3.1. PGL2(M) versus PSL2(R). A more familiar setup is that of automorphic 
functions on PSL 2 (R). We describe how it is connected to our setup. Let G\ = 
PSL 2 (R) and T\ C G\ be a (co-compact) lattice. We consider the quotient space 
X\ =T\\G\ ~ SY which is isomorphic to the tangent sphere bundle over a (com- 
pact) Riemann surface Y . The spectral decomposition of the Casimir operator on 
X\ (or of the Laplacian A on Y) leads to the decomposition L 2 (X\) = ®i(Ni,<Ji) 
into irreducible unitary representations of PSL 2 (R). Now consider Ti as a lattice 
in G = PGL 2 (R) (or in PSL^(R)). Formally, the space X = T\\G = X + U A_ 
consists of two connected components. It might happen however that there exists 
a lattice V d G such that Ti C T and rnPGLJ(I) 7^ 0. The best known example 
is Ti = PSL 2 (Z) and T = PGL2(Z). A geometric example could be constructed 
from a bordered Riemann surfaces Y with a totally geodesic boundary, by taking 
the connected sum of Y and Y along the geodesic boundary. 

We want to connect representations of G in L 2 (X) and those of G\ in L 2 (X\). 

First let us consider the situation where r = Y\. Let S = 1 j. We have 
X = X + U A_ with X + = Ti \ Gi, and A_ ~ X\ under the multiplication by 5, 
i.e., x i-> x5. 

Let (N, a) C L 2 (X\ ) be an automorphic representation of G\ which as an 
abstract representation is induced (i.e., of principal or complimentary series). 
Consider the space N' C L 2 (A_) of functions obtained by the action of S (i.e., 
N' = {/' I f'(xS) = f(x) for some / € N}). We have then N®N' = L ®L U 
where L £ , e = {0, 1} are (even and odd) irreducible representations of G which as 
representations of Gi are isomorphic. Under the restriction map L 2 (X) — > L 2 (Xi), 
these give the same space of eigenfunctions of the Casimir on X\ . 

For an automorphic representation (N, a) of G\ of discrete series, one first con- 
siders the complex conjugate representation (N,a) (i.e., the space of functions /, 
/ e N with the natural action of Gi), and then the action of 6 (interchanging 
connected components). In that case, we obtain two copies of isomorphic repre- 
sentations of G, and each of these split into two irreducible components under the 
restriction to G\ (holomorphic and anti-holomorphic discrete series representations 
on different components of X). 
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Taking restrictions of functions on X to subsets X± ~ X\, we obtain two 
isomorphisms i± : L 2 (X±) ~ i 2 (Xi) as representations of G\. Note that i± arc 
not coming from the representation-theoretic restriction of representations of G to 
representations of G\. 

The case of X consisting of one connected component is more interesting. Ge- 
ometrically this means that the oriented hyperbolic Riemann surface Y is iso- 
metric (as an oriented Riemann surface) to itself after reversing the orientation. 
We have X\ = X, and 5 acts on it on the right. Let us assume for simplicity 
that 6 normalizes Ti. Consider an automorphic irreducible unitary representation 
(L, ir) C L 2 (X) of G. If 7r = 7Ta, £ is induced, then we can take its restriction to 
the representation of G\ which is again irreducible. This leads to even and odd 
eigenfunctions (with respect to S) of C on X depending on the parity of e. Note 
that we can not reverse the direction by starting with some automorphic repre- 
sentation of Gi and extending action of G to it (if there are multiplicities, one has 
to split the space into the sum of <5-even and <5-odd subspaces). 

For an automorphic discrete series representation of G, the restriction to G\ 
leads to two irreducible automorphic components which are interchanged by 5. 

2.4. Geodesic circles. We start with the geometric origin of the spherical Fourier 
coefficients. We fix the "standard" maximal compact subgroup K — PO(2) C G, 
and the identification G / K — >• H, g h-> g ■ i, with the hyperbolic plane. Let 
T C PSL2(M) be as before. We denote by p-u :H^T\H^Y the corresponding 
projection. It is compatible with the distance function d(-, •) on Y and H. Let 
y e Y be a point, and let R y > be the injectivity radius of Y at y. For any 
r < R y , we define the geodesic circle of radius r centered at y to be the set 
cr(r,y) — {y' € Y\d(y',y) = r}. Since the map pu is a local isometry, we have 
that Puipui?-, z)) = a(r, y) for any z e H such that pn( z ) — where a-n(r, z) is 
the corresponding geodesic circle in H (all geodesic circles in H are the Euclidian 
circles in C, though with a different center from z). We associate to any such circle 
on Y an orbit of a compact subgroup on X. Namely, let Kq = PSO(2) C K be 
the connected component of K . Any geodesic circle on T-L is of the form a-u(r, z) = 
hKog-i with h, g € G such that h-i = z and hg-i € o~u(r, z) (i.e., an /i-translation 
of a standard geodesic circle centered at i € % and passing through g ■ i € H) . Note 
that the radius of the circle is given by the distance d(i,g ■ i) and hence g Kq 
for a nontrivial circle. Given the geodesic circle o~(r,y) C Y, we consider a circle 
a-u{r,z) C W projecting onto cr(r,y) and the corresponding elements g, h e G. 
We denote by K a = g~ l K$g the corresponding compact subgroup and consider 
its orbit JC a = hg ■ K a C X. Clearly we have p{JC a ) — a. We endow the orbit 
JCu with the unique ^-invariant measure d[i}c a of the total mass one (from the 
geometric point of view a more natural measure would be the length of a). 

We note that in what follows the restriction r < R y is not essential. From now 
on, we assume that K, C X is an orbit of a maximal connected compact subgroup 
Kq C G. The group K' is conjugate to PSO(2), but does not coincide with it. 
We denote by K' the corresponding maximal compact subgroup in G. Let JC C X 
be an orbit of Kq. Since we have assumed that X consists of two components 
X = X + U X- , we also have JC ■ K' = tC + U /C_ , JC± C X± and each of these sets is 
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an orbit of K' . We will assume that JC = K.+ C X+. The restriction r < R y simply 
means that the projection p{JC) C Y is a smooth non-self intersecting curve on Y . 
We remark that polar geodesic coordinates (r, 9) centered at a point z 6 % could 
be obtained from the Cartan XAiT-decomposition of PSL2(1R) (see [He] where 
analysis on H is discussed in the context of the representation theory) . 

2.5. Generalized periods and equivariant functionals. We fix a point o 6 K,. 
To a character x : K' — > C x , we associate a function x.(ofc') = x(k'), k' <E K' on 
the orbit 1C, and the corresponding functional on C°°(X) given by the generalized 
period 

(2-2) / /(*)*• 

for any / € C°°(X). The functional d^"^ is %-equivariant with respect to the right 
action of K' Q on functions on X: d^(R(k')f) = x{k')d a ^{f) for any k' € X£ 
(here R is the right action of G on the space of functions on X). For a given 
orbit /C and a choice of a generator xi of the cyclic group 1£q ~ Z of characters 
of the compact group -K^, we will use the shorthand notation d^ = d°^ K , where 
Xn — Xi - The functions (Xn). form an orthonormal basis for the space L 2 (K, dfx/c)- 

Let v : V — > C°°(X) be an irreducible automorphic representation. When it 
does not lead to confusion, we denote by the same letter the functional d° u ^ = 
^xK v on t ne space V induced by the functional defined above on the space 
C°°(X). Hence we obtain an element in the space Hom^(y,x). We next use the 
above mentioned multiplicity one property, i.e., the fact that the pair (G,K' Q ) is a 
Gelfand pair (i.e., that dim Horn^' ( V, x) < !)• 

Let V ~ Va, £ be a representation of the principal series. We have then the 
multiplicity one property dimHom^(VA, e , x) ^ 1 f° r am/ character x of Kq (i.e., 
the space of .K^-types is at most one dimensional for a maximal connected compact 
subgroup of G). In fact, dimHom^(VA,e, Xn) = 1 if and only if n is even. 

To construct a model %-equivariant functional on V\. £ , we consider the circle 
model V\ t£ — C^(S 1 ) in the space of even functions on S 1 and the standard 
vectors (exponents) e n = exp(in9) £ C°°(S' 1 ) which form a basis of K^-types for 
the standard connected maximal compact subgroup Kq = PSO(2). For any n 
such that dhnHomK (V\. £ , Xn) — lj the vector = ttx(<7 )e n defines a non-zero 
(Xn, i^o)-equivariant functional on Va j£ by the formula 

(2-3) C , W=Cl(») = < B '0 ■ 

We call such a functional the model x„-cquivariant functional on the representation 
V ~ V x , e . 

The uniqueness principle implies that there exists a constant a n = o, Xn xi l/ ) G C 
such that 

(2.4) C'M^n-rt), 

for any it £ V. Here we suppressed the dependence on the orbit and the auto- 
morphic representation since these are fixed in our discussion (nevertheless this 
dependence is central in other applications). 
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2.5.1. Invariant functional. Since the i^g-invariant functional will play a double 
role in our construction, we introduce another notation for it. Note that such a 
functional automatically vanishes on a representation of discrete series. 

Let xo = 1 be the trivial character of K . We denote by d T ^ £ {u) = {u,e' )v Tie , 
u G V TtS , the corresponding model functional. Let i/j : V n>ei — > C°°(X) be an 
irreducible automorphic representation of class one. We have as before 

(2.5) d? ttt (u) = d^ Ut {u) = f Ui(u)(k)xo.(k)dHK = a(i)d Ti , ei (u) , 

JK 

for any u g V Ti . £il and a constant a(i) = a/c(vi) e C. 

We want to compare the coefficients a(i) with more familiar quantities. Let /C = 
xq ■ Kq C X and let <p' T . — fi(e' ) be the automorphic function which corresponds 
to a i^Q-invariant vector e' Q € V n , Si of norm one. From the definition of d T _ ei it 
follows that 

(2.6) a{i) = <t>' ri {x ) . 

Hence, under the normalization that we choose, the coefficients a(i) are equal to 
values at a point xq for Maass forms on the Riemann surface Y' = T \ G/g~ 1 Kg. 

Remark. A trivial, but important, remark is that on the discrete series represen- 
tations any i^Q-invariant functional is identically zero, unlike a T-invariant func- 
tional for the diagonal subgroup T of G. This greatly simplifies the technicalities 
in what follows; and in fact, this is the reason that we treat only geodesic circles 
and leave closed geodesies aside. For the later, we need to discuss discrete series 
representations at length, and the resulting computations are more involved. 

2.6. First Gelfand pair: abelian spectral decomposition on a. Let v : 

V\. £ — > L 2 {X) be an automorphic representation of principal series. Consider the 
Hcrmitian form on V\ tS arising from the restriction to the orbit JC of Kq (i.e., to a 
geodesic circle a) 

(2.7) h\{4> u ) = j \<f) u \ J 2 d/C, 

JK 

for 4> u — v{u), u £ V\ t£ . Obviously the Hcrmitian form h\ could be decomposed 
with respect to the action of K' . We have 

(2.8) h x (ct> u ) = £ |dr 4 (^)| 2 - £ \a k \ 2 \dr d {u)\ 2 . 

k k 

This decomposition could be used in order to evaluate the integrated form 

(2.9) H\ = j k-h x dk . 

Jk 

This is the Hermitian form which appears in Theorem 11.41 restricted to an irre- 
ducible component inside of the space Ua ■ 

The form H\ is Ko-'mvaiia.Tit, and hence is determined by its values on the 
(essentially unique) orthogonal (with respect to H\) basis of V\ j£ consisting of K$- 
equivariant vectors {e„} ne 2Z described in section l2~2l Hence we need to estimate 
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the quantities H\((f> en ) = h\{4> en ), where we denote 4> en — v\{e n ). Substituting 
(|2.8p into the expansion (|2.6[) . we see that 

(2.10) H x (4> en ) = h x {<p en ) = ]T M 2 K lod (e n )\ 2 , 

feG2Z 

for any n 6 2Z. 

In order to prove bound (|1.5[) . we need to show that H\(cj) en ) <C \n\ s for any 
S > and fixed A. 

2.7. Theorem. For any S > 0, £/iere exists a constant A\ 7 g > smc/i </ia< 

(2-11) ]>> fe | 2 |dr d (e„)| 2 <A M M 5 , 

fe 

/or any n € 2Z. 

This implies bound (|1.5p and proves the main theorem of the paper. Below 
we present the proof for the principal series representations, and discuss in [X] 
amendments needed to be made for the discrete series. 

Remarks. 1. The coefficients ak are well known in the theory of automorphic 
functions. These are called the spherical Fourier coefficients of Maass forms and 
were introduced by H. Peterson. Quantities |afc| 2 are related to special values of L- 
functions (for special geodesic circles on Riemann surfaces of arithmetic origin, and 
for what is called Hcckc-Maass forms) , and are of utmost importance in Number 
Theory. In particular, it is generally believed that they satisfy the bound \ak\ <C 
\k\ e for any e > (the Lindcloff conjecture). We stress that coefficients ak depend 
on fj, v and T. 

2. It is relatively easy to prove the mean- value bound 

(2.12) $]k| 2 <AT, 

fc|<T 

for any T > |A|, and some universal constant A depending on Y only (this follows 
from |BR2j . and was spelled out explicitly in |Re2) ) . 

2.8. Oscillating integrals. We treat principal series representations first and 
discuss discrete series in [3] The coefficients d 1 k nod (e n ) should be viewed as a func- 
tion of two variables k and n (in fact, there is also a dependence on the parameter 
of representation A which we suppress as it is fixed in our discussion). Recall that 
we defined these by the matrix coefficient (|2.3[) 

(2.13) dr d (e n ) = (e n ,n x (g~ 1 )e k ) . 

Here ej £ V\ :E is the norm one j-ifo-type given by the function ej(9) = e ye in the 
circle model Vx,e — G^^S 1 ). Such a matrix coefficient clearly could be computed 
via an oscillating integral in one variable. Taking into account the action of Kq in 
the space V\, e , we have 

(2.14) d k nod (e n )= [ \ g '(8)\-h^ h ^ 9 '^e lk9{e) - me d6 . 

Js 1 

Note that this formula does nor depend on the parity e of V\ y£ . This is not 
surprising as e describes the action of the element 5. 
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Integrals (|2.14l) are common in the theory of special functions, and have been 
well-studied in Classical Analysis a la Whittaker and Watson. In fact, these are 
well known Legendre functions. The asymptotic of such an integral is controlled 
by the stationary phase method. We are interested in the behavior of d™ od (e„) 
for large parameters k and n, for A and g which are fixed. Hence we denote by 
A{6) = |. 9 '((9)|-5e5 Aln l9'( e )l the amplitude and by S k . n (9) = kg(9) - nO the phase 
in the oscillating integral (|2.14p . By computing the relevant phase function, one 
can easily see that the function d™ od (e n ) is well modeled by the classical Airy 
function A(a;) (see |Hel Theorem 7.7.18]; we note that careful analysis of such a 
reduction is done in |BR3j ). 

Let M = max eeS i \g'{9)\ > 1. We have M' 1 < g'(9) < M as g is a linear- 
fractional map on S 1 . We have three types of behavior for the integral (|2.14l) . 

(1) For k in the range lAM~ 1 n < k < 0.9Mn, the phase function has two 
critical points of Morse type; in this case we can estimate the integral 
using the stationary phase method. In this range the integral is of order 
of \k\~i. 

(2) For k close to M~ 1 n or to Mn, critical points of the phase collide to a 
cubic critical point. In order to get uniform bounds in this region, we can 
use properties of the Airy function. 

(3) For k in ranges k < 0.9M _1 n or k > lAMn, there are no critical points 
of the phase, and the value of the integral is negligible (in n). This region 
could be discarded in further analysis. 

We treat the case k x Mn, and the complimentary case k x M~ lr n is identical. 
Namely, we assume that n > and k > n (the treatment for the range k < n is 
identical). Let us introduce the notation 6 — S(k, n) — Mnjk. 

The above analysis (made rigorously in |Re3) and |BR3j ) implies in particular 
that for any fixed sufficiently small e > 0, in the region 1 + e > <5 > 1 — e (i.e., 
k x Mn), we have 

(2.15) \dr d (e n )-A(k,S)\ <C\k\~i . 

Here A(k, 5) = \k\~%A(k*(5—i)) = \k\~TsA(k~i(Mn— k)), where Ais the classical 
Airy function (see [He]). We only need to know that the Airy function A(x) is a 
smooth function, bounded at 0, rapidly decaying for x > and uniformly bounded 
by |x| -1 / 4 at infinity. This well-known asymptotic of the Airy function imply that 
(roughly) we have three types of behavior: 

Regular: If M - 1 < (Mn - k)/k < 1 - e, then K" od (e„)| < |n|~i 
Resonance: If 1 - s < (Mn - k)/k < 1 + e, then \df od (e n )\ < \k\-?(l + 

|jb| 2 / 3 |Mn-Jb|)-V* 
Cutoff: If (Mn -k)/k>l + e, then \d™ od (e n )\ < \k\~ N for any N > 0. 

Splitting the summation in (|2.11l) according to these regions of k, we see that 
the mean- value bound (|2.12p allows one to treat the regular part of the summation 
in (I2.11[) , i.e., k such that S is not close to 1. This region of summation in k is of 
length comparable to n due to the cutoff, and the value of \d™ od (e n )\ 2 is of order of 
|ri| _1 . Hence we see that the total contribution from this region of k is uniformly 
bounded as n — » oo. 
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We are left with the "short" sums of the form Y^\s(k n)-i|<T- 1 \ a k\ 2 f° r n an d 
T — > oo, T -C n. Denoting T = Mn, we consider a sum X)||fc|-T|<|T|^ \ a k\ 2 f° r 
7 < 1. The size of the transition region of the Airy function implies that we only 
need to consider the range 1 > 7 > 2/3 (in fact, we do not know how to deal 
with shorter sums!). It turns out that these sums could be bounded effectively 
with the help of the second Gelfand pair (i.e., using the multiplicity one statement 
dim Home (V\ <S> V2 <8> V3, C) < 1). This is the more "tricky" part of the proof and 
is done in [Re3j . We explain it in the next section. 

3. Triple product 
Here we explain how to obtain the bound 

(3.1) J2 M 2 ^ B ^ +e . 

|fc|-T|<T§ 

This is the content of Theorem 1.5 from [Re3j . We recall the corresponding setup. 

3.1. Second Gelfand pair: triple product spectral decomposition. Recall 
that we consider the Hermitian form given by h\(4>) — L- | </> u |k;| 2 cWC for u £ V\ tS . 
We now switch to the corresponding linear functional on E — V\. e <E) V-\ tS which 
we denote by the same letter. We have h\(w) = J K r/s L {vE{w))dlC for w € E. We 
now consider the spectral expansion for the form h\ coming from triple products. 

3.1.1. Spectral theory on X. Let v : V — > C°°(X) be an irreducible automorphic 
representation as before and ve = v®v: E = V®V— > C°°(X x X) the 
corresponding realization of E. We assume that the space X is compact, and 
hence we have the discrete sum decomposition (I2.1[) 

(3.2) L 2 (X) = (®i(L h Ui)) © (® K (L K , v K )) 

into irreducible unitary representations of G. Here Vi : Li —¥ L 2 (X) are unitary 
representations of class one (i.e., those which correspond to Maass forms on Y 
with the eigenvalue (1 — rf )/4) and L K are representations of discrete series (i.e., 
those which correspond to holomorphic forms on Y). We denote by V; C 
the corresponding spaces of smooth vectors and by pr L . : L 2 (X) — > Vi(Li) the 
corresponding orthogonal projections (note that pr L . : C co {X) — > Vi). 

We consider triple products of eigenfunctions. Let r\ : C co (X x X) — > C°°{X) 
be the map induced by the imbedding A:I->IxI. Let Vi : V Tu£i — > C°°(X) 
be an irreducible automorphic representation. Composing r& with the projection 
pr Li : C°°(X) -> Vi(V Ti>ei ), we obtain AG-invariant map T? ut : E -> V TitSi and 
the corresponding automorphic trilinear functional lf ut on E ® V*. e . defined by 
I1 ut (v ® u ® t) = (r&(vE{u (g> v)), Vi(t)) (here we identified V Tiy£i with the smooth 
part of V*. e ~ V— TitSi ). Such a functional is clearly G-invariant, and hence we can 
invoke the uniqueness principle for trilinear functionals (see section f2 . 1 [) . 
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3.1.2. Triple product spectral expansion. Recall that in section r2.5.1l wc denoted 
by d^ ut : Vt — > C the i-Cp-invariant functional coming from the integration along 
the orbit K, C X (sec (|2.5p ). The expansion (|3.2I) implies that 

(3.3) h x ( W ) = Y / dr t (Tr t M) , 

i 

for any w £ E. We now write d^ ut (T"" f (w)) as an integral transform in the circle 
model of representations. 

3.1.3. Model triple product. Let V\^ e and V T>e > be two irreducible unitary repre- 
sentations of principal series (and we denote by E the smooth part of the tensor 
product V\,£<8> V-\_ e ). We consider circle models of these representations where the 
group K' acts by the standard rotations of S 1 . We construct a (non-zero) explicit 
functional l^^y , £ Home (E <g) V- T ^> , C) which we call a model functional. Let 
us denote by sgn(9,9',9") = sgn ((9 - 9') (9 - 6") (6' - 9")) the function taking- 
values ±1 on (S 1 ) 3 \ {points with at least 2 coordinates equal}, and changes the 
sign at points removed. There are exactly two open orbits for the diagonal action 
of PGL^(M) on (5* 1 ) 3 , and the above function is an invariant of these. Moreover 
it is antisymmetric with respect to the action of the element 6. It is easy to see 
(as in [BR2] ) that in the circle model of class one representations the kernel of 
tfstsv i i s gi ven by the following function in three variables 6. 9' , 9" £ S 1 

(3.4) X A ,_ A , r (M',#") = (sgn(9,9',9")y' '■ 
| sm(9 -9')\^\sm(9-9")\ \ sm(9' - 9")\ . 

(The factor should be equal to (sgn(9,9' ,9")) e+e+e , but it gives the same func- 
tion.) This function also defines the kernel of the map via 
the relation 



(T T (w),v) v ,=t^/ w{9 : 9')v(9')K X - XT (9,9',9")d9d9'd9'' 
/(si)3 



(2tt) 3 



3.1.4. Integral transform. The model functional d T is given as the scalar product 
with a norm one /^-invariant vector e' (9") = 1 (i.e., d T {v) ~ (f,eo)v re , as m 
section [2.5. ip . Hence we have 

(3.5) d T (T T (w)) — (T T (w),e' )v T — 

w{9, 9')K X -\A e i 9 'i 9")e' {9")d9d9' d9" , 



(2tt 



for any w £ C°°(S 1 x S 1 ). For what follows, it will be enough to assume that the 
vector w £ E is A/l^-invariant. Such a vector w can be described by a function 
of one variable; namely, w{9,9') = u{c) for u £ C 00 ^ 1 ) and c = (9 - 9')/2. We 
have then w(n, —n) — u(n) = ^ J sl u{c)e~ mc dc - the Fourier transform of u. 
We introduce a new kernel (note that e' Q (9") = 1 in (|3.5p ) 

(3.6) k T (c) = fc A , T (Vi) = i / #a,-A,t(0, 0', 9")d9" 

2?r J S i 
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and the corresponding integral transform 

(3.7) uHt) = «» (r) = — ^ / u(c)k T (c)dc , 

(lliy Jgl 

suppressing the dependence on A as we have fixed the representation V\ t€ . The 
transform is clearly defined for any smooth function u € C°°(S 1 ), at least for 
t G iM.. In fact, it could be defined for all r S C, by means of analytic continuation. 
We will discuss this in \K\ where we deal with discrete series. 

Note that k T is the average of the kernel K\.-\ :T with respect to the action of 
AK' , or, in other terms, k T — T*(e' Q ) € E* is the pullback of the i^-invariant 
vector e' € V T ,e' under the map T*. 

3.2. Gelfand formation: spectral identity. We invoke now the uniqueness 
principle dim Home (E ® V^ T _ £ ',C) < 1. This implies that the two invariant 
functionals, lf ut and l%g d v _, , we constructed are proportional. Hence there are 
constants — fi{vi) <E C such that 

(3.8) T z aut = (3{i) ■ T T( . 

We also recall that the automorphic if -invariant functional on V Ttet satisfies the 
relation d^ ut — a(i) ■ d Ti as in (|2.5.ip . Denoting by — a(i)(3(i) the product of 
these constants, we rewrite the spectral expansion (|3.3[) in the form 

(3.9) h x (w u ) = Y J dr t {Tr\w)) = ^ 7 «rfn(T n K)) 

Ti 

for any Ai^Q-invariant vector w u £ E corresponding to the function u € C°°(S 1 ). 
This is the spectral expansion corresponding to the triple product. Taking this 
together with the abelian spectral expansion (|2.8I) , we obtain the spectral identity 
corresponding to two Gelfand pairs (we call such a pair a Gelfand formation, and 
the corresponding identity the period identity): 

(3.10) \ a n\ 2 Hn) = u(0) + 7(0 ' ■ 

Here on the right we singled out the contribution from the trivial representation 
(i.e., r = 1), and wrote it in the form u(0) = (Vol(/C)/ Vol(X)2) • u(0) under our 
normalization of measures Vo\(X) — 1 and Vol(/C) = 1. 

By choosing the appropriate test function u, and estimating the right hand side 
in (|3.10[) . we will obtain the bound (|3.1I) . 

3.3. Bounds for spherical Fourier coefficients. We are interested in getting 
a bound for the coefficients a n . The idea of the proof of the crucial bound (|3.1j) is 
to find a test vector w £ V <E) V, i.e., a function w € C 00 (S' 1 x S 1 ), such that when 
substituted in the period identity (|3.10l) it will produce a weight w which is not 
too small for a given n, \n\ — > oo. We then have to estimate the spectral density of 
such a vector, i.e., the transform wK One might be tempted to take w such that 
w is essentially a delta function (i.e., picks up just a few coefficient a n in (|3.10|1 ). 
However, for such a vector we have no means to estimate the right hand side of 
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the formula (|3.10p because w$ is spread over a long interval of the spectrum. The 
solution to this problem is well known in harmonic analysis. One takes a function 
which produces a weighted sum of the coefficients |afe| 2 for k in a certain range 
depending on n and such that its transform vfi is spread over a shorter interval. 
For such test vectors w, we give an essentially sharp bound for the value of the 
diagonal period cIaic(w) — § K r&(vE(w))dlC. 

3.3.1. Proof of bound (|3.1[) . We start with the period identity p.lOp and construct 
an appropriate AA'p-invariant vector w G E, i.e., a function u € C°°(S 1 ) such that 
w(8,8') — u((9 — 6')/2). We have the following elementary technical lemma. 

3.4. Lemma. For any integers N > T > 1, there exists a smooth function u^.t G 
C^iS 1 ) such that 

(1) K,t(0)I < aT, 

(2) u N ^ T {k) > for all k, 

(3) UN,r{k) > 1 for all k satisfying \k — N\ < T , 

(4) |< T (r)| < aT\N\-?(l + \r\)-i + oT(l + |r|- 5 / 2 ) for \r\ < N/T, 

(5) |< T (r)| < oT(l + |r|)- 5 / 2 for \r\ > N/T, 
for some fixed constant a > independent of N and T. 

The proof of this Lemma is given in Appendix A of |Re3j . One constructs 
the corresponding function un.t by considering a function of the type un,t{c) = 
rp e -zNc . ^ (t c ^ w jth a fixed smooth function -0 G C°°(S 1 ) of a support in a 
small interval containing 1 e S 1 (here * denotes the convolution in C°°(S 1 )). Such 
a function obviously satisfies conditions (1) — (3) and the verification of (4) — (5) 
is reduced to a routine application of the stationary phase method (similar to our 
computations in |BR3j ). □ 

We return to the proof of bound (|3.ip . In the proof we will use two bounds for 
the coefficients a(i) and j3(i). Namely, it was shown in [BR2 that 

(3.H) £ \m\ 2 < ^ , 

A<\n\<2A 

for any A > 1 and some explicit a > 0. 

The second bound that we will need is the bound 

(3.12) J2 Hi)\ 2 <bA 2 , 

A<\n\<2A 

valid for any A > 1 and some b. In disguise this is the classical bound of L. 
Hormander [H] for the average value at a point for eigenfunctions of the Laplace- 
Beltrami operator on a compact Riemannian manifold (e.g., A on Y). This follows 
from the normalization |/3(«)| 2 = \4>' T . (a;o)| 2 we have chosen in (|2.6p for A~Q-invariant 
eigenfunctions. In fact, the bound (|3.12p is standard in the theory of the Selberg 
trace formula (see [Iwj ) . and also can be easily deduced from considerations of 
[BR3] , 

We plug a test function satisfying conditions (1) — (5) of Lemma [3.41 into the 
identity (|3.10p . Using the Cauchy-Schwartz inequality and taking into account 
bounds (|3.1ip and (|3 . 1 2|) . we obtain 
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\ a k\ 2 < J2 \ak\ 2 u N ,T(k) = u N . T (0) + a(0/3(*) u 5vT( T ») 

\k-N\<T k n^l 

<aT+ J2 aT\N\-h(l + \Ti\)-i\a(i)l3(i)\+ £ a<r(l+\n\)-*'*\a(i)0(i)\ 

\Ti\<N/T Tj^l 

<aT + aT|7V|^ £ (1 + | T< |)-* (|a(i)| 2 + |/3»| 2 ) + 
|n|<JV/T 

aT X) (l + |Ti|)- 5/2 (l«(*)| 2 + |0(*)| 2 ) <aT + CT|7V|-^(iV/T) 3 / 2 + £ + OT = 

r,#l 

c , T + c , T -i- £ | iV |l + e : 

for any e > and some constants c', C, D > 0. 

Setting T = N 2 / 3 , we obtain Y,\k-N\<W/* K| 2 < A £ N 2 / 3+E for any e > 0. 

This finishes the proof of the bound (|3.1j) . and with it the proof of Theorem 
12.71 for representations of principal series. In [A] we explain the case of discrete 
series. □ 



Appendix A. Representations of the discrete series 

The aim of this appendix is to indicate changes needed in order to treat discrete 
series representations. The main difference in treatment of the discrete series is 
the lack of their convenient realization. The problem is computational and not 
conceptual. We show how to reduce necessary computations to the case of induced 
representations. With these changes, the treatment is identical to the case of 
induced representations we discussed in the main body of the paper. The main 
reason why we can make this reduction comparatively easy is that we assume that 
the representation is fixed. Similar computations that take into account the weight 
of the representation would be more complicated. 

A.l. Representations and their realizations. Let k > 2 be an even integer, 
and (Dk,TTk) be the corresponding discrete series representation of PGL^M). In 
particular, for m £ 2Z, the space of Ko-tyj>es of weight m is non-zero (and in this 
case is one-dimensional) if and only if |m| > fc. This defines -Kk uniquely. Under 
the restriction to PSL^R), the representation Wk splits into two representations 
{p^,wt) of "holomorphic" and "anti-holomorphic" discrete series of PSL^R), 
and the element 8 interchanges these. 

We consider two realizations of discrete series as subrepresentations and as quo- 
tients of induced representations. Consider the space Hk-2 of smooth even homo- 
geneous functions on R 2 \0 of the homogeneous degree k — 2 (i.e., f(tv) = t k ~ 2 f(y) 
for any t £ R x and / d £ I 2 ). We have the natural action of GL^R) given 
by ^k(g)f(v) = f{g~ 1 v) ■ det(g)( fc ~ 2 )/ 2 , which is trivial on the center and hence 
defines a representation of PGL^R). There exists the unique non-trivial invariant 
subspace 14_2 C Hk-2- The space 14_2 is finite-dimensional, dimVfe_2 = k — 1, 
and is generated by the homogeneous polynomials of degree k — 2. The quotient 
space "Hfc_2/Vfc_2 is isomorphic to the space of smooth vectors of the discrete series 
representation nk- Note that Vk—2 consists of vectors with ifo-type in the range 
| raj < k. 

We also consider the dual situation. Let T-L-k be the space of smooth even 
homogeneous functions on R 2 \ of the homogeneous degree —k. There is a 
natural PGL2(R)-invariant pairing ( , ) : Hk-2 ®H—k — > C. Hence H-k is the 
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smooth dual of Hk-2, and vice versa. There exists the unique non-trivial invariant 
subspace in H-k, and it is isomorphic to Dk- The quotient H-k/Dk is isomorphic 
to the finite-dimensional representation V k -2- 

Consider the restriction of smooth functions in Hk-2 to functions on the circle 
5* 1 C K 2 \ 0. There exists a unique (up to a multiple) Hermitian G-invariant form 

n— — k n— oo 

onHk-2 ^ C^iS 1 ) givenby F k (f,g) = E J-(k,n)-a n b n + l+(k,n)-a n b n , 

n=-oc n=k 

where a n and b n are the Fourier coefficients of / and g correspondingly, and the 
coefficients ~f±(k,ri) are given by 7±(fc,n) = T(k ± n)/T(±n — k). (Note that 
Fk is degenerate on Hk-2-) The Hermitian form Fk corresponds to the invariant 
Hermitian product on 7Tfc after passing to the quotient (in particular, it is positive 
definite on the quotient). We denote it ( , )jj k . Note that asymptotically the 
coefficients r )±(k,n) grow as |n| 2fe+1 , and hence the corresponding norm on irk 
resembles the Sobolev norm on S 1 . We will not use the invariant norm, and 
instead use the natural linear pairing on Dk ® F>* k and its relation to the natural 
pairing on Hk-2 ® H-k- 

Consider the representation Dk and its (smooth) dual D* k (of course Dk — D k , 
but we want to distinguish between these two copies). We want to find a way to 
compute the canonical pairing ( , )k ■ Dk <8> D k — > C in terms of the canonical 
pairing ( , ) : Hk-2 ® H-k — > C. Consider the imbedding ik : D k — > and 
the quotient map tfe : 'Hfe_2 — > -Dfc- Hence we also have id® ik ■ Hk-2 ® — )• 
Hk-2®H-k and qk®id : H k -2®H-k -> D k ®H- k - The image of the composition 
qk®idoid®ik is equal to the image of the imbedding id®i k ■ Dk®D k — > Dk®H-k- 
Hence for a pair of vectors u <S> v E Hk-2 ® such that w = ik(v'), v' G D^, we 
have (qk(u),v)k — (u, v'). This will be our way to compute the invariant Hermitian 
norm on Dk- The important difference, as opposed to the invariant Hermitian form 
( , ) o k on Dk, is that it is not possible to say what is the norm of a vector v £ Dk 
(for this, one uses the intertwining operator between H-k and Hk-2 which we want 
to avoid). It is however possible to say when the scalar product between a vector 
and a co-vector is one. In general this does not allow one to construct a projector 
to a particular vector. However in our situation we will need projectors onto pure 
-Ko-types, and since these are essentially unique the above construction is enough 
for our purposes. We note, however, that for a general pair u <E> v € Hk-2 <8> H-k 
not satisfying the condition v — ik(v') (i.e., v £ Im(ik) C H-k), the value of (u, v) 
might be far from a pairing of vectors in the discrete series (e.g., paring of high 
-Ko-types components of u, v). This is because components belonging to the finite 
dimensional representations might be dominant in the pairing ( , ). 

We now compare this to the automorphic picture. Let v : Dk — > L 2 (X) be an 
isometric automorphic realization of a discrete series representation. Consider the 
complex conjugate realization v : D* k — > L 2 (X). Integration along the diagonal 
AX X x X corresponds to the canonical automorphic pairing ( , ) k ut : Dk® 
D k — > C. We have ( , ) k ut — ( , )k since it is assumed that v is an isometry. 
While as we noted it does not allow us to compute norm of a vector in Dk, we can 
produce some functions on X corresponding to norm one vectors. For example, 
let e„ C Dk be a vector of i^o-type n, and e*_ n G D k the dual if -type vector (i.e., 
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(e n ,e*_ n )k = 1)- We have then v(e n ){x)v{e*_ n ){x) = \v(e n )(x)\ 2 as functions on 
X, where e n G Dk is a norm one i^o-type. 

A. 2. Invariant functionals. We now construct i^o-equivariant functionals and 
triple functionals for discrete series using the space Hk-2 <8) H-k- Let v ® v : 
Dk®D^ — > L 2 (X x X) be as before, and K. C X be a i^-orbit. For a character % : 
K' — > S 1 , we define the automorphic functional = f K v(u)(k)xn(k)dJC 7 

which is proportional to the model functional d n (u) = (u,e' n )D k , i-e., dj"* = 
a n ■ d n , where e' n € Dk is a norm one vector of i^-type n. On the other hand, 
we have the corresponding functional d n <g> d-„ : %fc_2 ® H-k —> C given by 
d„ <g> gL„(u ® u>) = (it, e' n ) ■ (w, e'_ n ), where ® e'_„ e H-k ® Hk-2 is any pair of 
corresponding ii^-types which are in duality (e' n ,e'_ n ) = 1. The main point here 
is that we can easily compute the value of d n ®d- n (u®w) since in the circle model 
H-k®Hk-2 — C co (S 1 xS 1 ) these are Fourier coefficients of functions u and w. For 
\n\ > k, we have ^y"*®^x-n = |an| 2- dn<£>d-n- Let ej € be a norm one j-th Kq- 
type (we recall that there are two different maximal compact connected subgroups 
K = PSO(2) and K' Q involved in the description of the set Ocl). As in (|23Uj) . 
we are interested in coefficients \d n (ej)\ 2 — \{ej, e' n ) u k \ 2 ■ For \j\ > k, a vector with 
the K -type equal to j automatically belongs to the subspace Dk C H-k, and we 
have \d n (ej)\ 2 — d n {ej) ■ d- n (e—j) = d n {e.j) ■ dL n (e_j-). Here e±j € Dk are norm 
one vectors and e±j € "H are any vectors in duality (e"j,e'_j) = 1. This implies 
that we have the integral representation for the quantity similar to the integral 
(|2.f4p we used for principal series. In particular, this function exhibits the same 
Airy type behavior. Hence as for the principal series, we need to estimate the 
coefficients a n . For this we construct model triple products as in the case of the 
principal series. 

Let ( V Tt£ i , TT T ,e' ) be a unitary representation of the principal series. The space 
Home (£>fc<g)-D£, V T . e i) is one-dimensional. We will work with the space of invariant 
trilincar functionals Homc(Dk®D' k '(&V- T , e ', C) instead. We first construct explic- 
itly a non-zero element lk, T ,e m the space Homc(^fc-2 <8> H-k ® V- T:€ >, C) (which 
is in fact also is one-dimensional), and then use it to define a non-zero element in 
the space Hom^D^ (g> Df. ® V- T ,e', C). What is more important, we will use lk, T ,e 
to carry out our computations in a similar way to the principal series. 

Consider the following function (compare to (|3.4[0 in three variables 9, 9' , 9" £ 

S 1 

K k - 2 ,-kA e ^'^") = {sgn(9,e',e")f ■ 

I sin(# - 9')\^\sin(9 - 9")\^' k+1 \ sm{9' - 9")\^+ k ~ 1 . 

Viewed as a kernel, it defines an invariant non-zero functional on the (smooth 
part of) the representation Hk-2 ® H-k ® V- T>e > ~ C^S" 1 x S 1 x S 1 ). Such a 
kernel should be understood in the regularized sense (e.g., analytically continued 
following |Glj ). However, since we are interested in r 6 iR, all the exponents in 
(|A.I[) are non integer. This implies that the regularized kernel coincides with the 
above function when integrated against any test function vanishing in a neighbor- 
hood of singularities of the kernel (lA.f I) . We denote the corresponding functional 
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by lk,T,e € HoniG(Hjt_2 S3 W-fc ® V1 T)£ <,C). Such a functional defines a non-zero 
functional lk,r,s G HomcC-Dfe ® ® y_ Tje <,C) by requiring that Zfc, T , e and Zfc, T , e 
induce the same functional on the space Hfc_2 <8> D% ® F_ T)E ' . We call Zfc iT , e the 
model functional for the discrete series. The difference with principal series clearly 
lies in the fact that we only can compute the auxiliary functional lk, T ,s- However, 
for k fixed, it turns out that necessary computations are identical to the ones we 
performed for the principal series. 

As in we are interested in fc fcjT {^f) = ± J gl K k ^ 2 ,-k,r(0,9' ,6")d9" . 

Moreover, we are interested in asymptotic of this function as |r| — ¥ oo and k is 
fixed. The analysis of such an integral is standard in the theory of the stationary 
phase method, and is similar (since k is fixed) to one which is done in IRe3] and 
[BR3 . This is because the contribution from any small enough neighborhood 
of the singularity of the above kernel is negligible when |r| — > oo due to the high 
oscillation of the kernel. The situation is of course different if also k — > oo, but this 
is not relevant to our present situation. Hence we obtain the asymptotic expansion 
of the trilincar invariant functional essentially identical to the one we obtained for 
principal series. This allows us to use the same calculations which we already 
explained before. In particular, we have to construct test vectors analogous to the 
ones obtained from functions um.t € C°°(S 1 ) in Lemma [3.41 We can take the 
same function as for the principal series, but while it gives a Aifo-invariant vector 
in the space T-Lk-2®'H-k such a vector does not belong to Hk-2 As a result, 

it is difficult to estimate the norm of such a vector. However it is easy to correct 
the function uat,t so that it will have vanishing Fourier coefficients un,t(ji) for 
\n\ < k (and hence will produce a vector in 1Lk~2 ® D* k ). Since we are interested 
in the asymptotic as N, T — > oo, this finite correction does not affect it. Here 
the fact that we consider fixed k is crucial. Without this restriction, it is more 
difficult to calculate the corresponding value of the trilinear invariant functional 
(an in fact the resulting asymptotic should be different). 
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